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Abstract. We calculate the effect of the Earth- Moon (EM) system on the free- 
fall motion of LISA test masses. We show that the periodic gravitational pulling 
of the EM system induces a resonance with fundamental frequency 1 yr~^ and 
a series of periodic perturbations with frequencies equal to integer harmonics of 
the synodic month (~ 3.92 X 10"'' Hz). We then evaluate the effects of these 
perturbations (up to the 6*'' harmonics) on the relative motions between each 
test masses couple, finding that they range between 3 mm and 10 pm for the 2"'* 
and G**" harmonic, respectively. 

If we take the LISA sensitivity curve, as extrapolated down to 10~® Hz in [T], we 
obtain that a few harmonics of the EM system can be detected in the Doppler 
data collected by the LISA space mission. This suggests that the EM system 
gravitational near field could provide an additional crosscheck to the calibration 
of LISA, as extended to such low frequencies. 
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1. Introduction 

LISA (Laser Interferometer Space Antenna) is a ten years long NASA-ESA space 
mission to detect gravitational waves in the frequency range 10~^ — 10"^ Hz [2]. It 
consists of three spacecrafts whose mutual distances are about L = 5 x 10^ km. The 
LISA constellation will orbit around the Sun following the same path of the Earth, 
(/)o = 20° behind [3]. 

The ideal configuration for LISA performances should be a rigid cquilatcr triangle 
[3|; however, the shape of the LISA constellation is subject to significant variations 
because of the gravitational interaction due to Sun and planets [S] . 
Much smaller perturbations can be induced by the presence of interplanetary dust [S] 
or dark matter in the solar system [7]. The perturbations due to each celestial body 
can be treated, at a first approximation, independently. The gravitational effects 
are quite different in intensity (orders of magnitude) and/or in behavior (stationary 
or time-dependent); the frequencies involved are, in general, not commensurable 
and so resonance effects are not observed. The only significative exception is the 
Earth perturbation [5], [8] and [9], which gives a resonance because of the 1:1 
commensurability between the Earth and LISA orbits. 

However, such a resonance will not have enough time to grow during the 10 years 
of LISA mission. In this paper we focus on the perturbations induced by the EM 
system on LISA at the frequency of the synodic month and its harmonics, which 
are much higher than 1 yr~^ and so, in a first order approximation, perturbative 
effects can be treated independently. We extend the approach of [5] to include the 
time dependent perturbation of the Moon. The plan of the paper is as follows. In 
Section[2]we shortly describe the perturbative approach used to study the EM effects. 
Section |3] is devoted to illustrate the EM system and the approximations we used to 
model its gravitational near field. In Section 2] we calculate the modulations of the 
distance between two LISA test masses due to EM system. In Section [5] we estimate 
the perturbations induced by Venus and Jupiter, compared to those due to the EM 
system calculated in the previous Sections. In Section [6] conclusions are drawn and 
the future research potential of method for further studies and application to LISA 
are given. 

2. Perturbative dynamics of the LISA test masses 

In order to simplify the notation, we use the Astronomical System of Units, for length 
(AU), mass (Mq) and time (days). However, the quantities that affect the relative 
motion of the LISA spacecrafts will be reconverted in the SI units. 
By means of a F77 code based on the inverse 15*''-order Runge-Kutta method |10| . 
we calculate the modulus of the force gradient between 2 LISA spacecrafts due to the 
main Solar System bodies, i.e. Sun, Venus, Earth (more precisely the EM system) 
and Jupiter (Figured]). The effects are different by orders of magnitude in amplitude 
and show also different frequencies. 

In fact, the Earth and Jupiter tidal effects are 10"* and 10^ times smaller than the Sun 
contribution, respectively; after 3 years. Sun, Earth and Jupiter cause arms length 
changes (the so-called "arm breathing") of ~ 1.2x10^ km, 4.8 x 10^ km and 4x10'^ km, 
respectively [4], [5]. 

The Venus contribution oscillates in time by two orders of magnitude and it is 
comparable with Earth effects for a short time interval every 584 d [11]. However, 
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as we will show in Section [51 this perturbation is negligible with respect to the Earth 
one. 



0.0001 

< 
3 



~'. 10-8 
m 
3 
C 
ta 

> 10-' 
x: 

to 

M 10-8 
C 

^ io-» 

<] 

10-10 

2 4 6 a 10 

time (a012+yrs) 

Figure 1: Gravity gradients between two spacecraft due to EM system (solid line), 
Venus (dotted line), Jupiter (short dashed line) and Sun (long dashed lines). 

Wc therefore conclude that the celestial bodies affect the LISA arms length by at most 
a few percent. As a consequence, at any time the distance of each satellite from the 
LISA barycenter is po = L/\/3 ~ 1.9 x 10^^ AU, and the distance between the Earth 
and the LISA barycenter is rg = 2 sin((/)o/2) ~ 0.347 AU, within a relative fluctuation 
of few percent. 

2.1. Hill- Clohessy- Wiltshire (HCW) reference frame 

To study the EM system effect, we made use of the Hill- Clohessy- Wiltshire (HCW) 
reference frame {x,y,z} [12] defined as follows (see Figure [2]): 

(i) the origin O' describes a circular orbit on the ecliptic plane at 1 AU from the 
Solar System Barycenter O [T3] : 

(ii) the xy plane coincides with the ecliptic; 

(iii) the z-axis is perpendicular to the ecliptic plane and parallel to the Solar System 
total angular momenturrQ; 

(iv) the i-axis is tangent to the orbit and is antiparallel to the origin O' velocity 
vector; 

(v) the y-axis is directed radially outward. 

In this frame the coordinates of the EM system barycenter are {xg,yg,Zg) = 
Ro{— sin 00, cos 0o — Ij 0). Of course, the LISA spacecrafts can be considered as three 




tt The Solar System total angular momentum is perpendicular to a piano slightly inclined with respect 
to the ecliptic one. For our scopes, such inclination is negligible. 
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Figure 2: Hill-Cloliessy- Wiltshire frame: Rq is the radius of the orbit, p and R are 
the the distances of P from the origin of the rotating and fixed frames, respectively. 
{X, Y, Z} is the associated reference frame (see details in Section [3]). The EM system 
is also represented. The figure is not to scale. 



proof masses. Their equations of motion in the HCW frame read 



.2 1^ 
2woy - ujqX = 



y + 2uJox-uj^iy + Ro)^--^{RQ+y) (1) 
.. /f_ 

w here ^ ^ GMq = uj^R^, ujq = 27r/365.257d-\ Rq = I AU , and R = 
-^/x^ + (y + i?o)^ + z^- Being the breathing length Ai ^ L, each satellite is located 
at any time at a distance p{t) = x{t)^ + y{t)^ + z{t)^ ~ po- Since po <C i?0 7 we 
expand the acceleration due to the Sun in terms of x/ Ro,y / Rq, z/ Rq. Retaining the 
first term of the series, we have the so-called Hill- Clohessy- Wiltshire equations [12] of 
relative motion 

X — 2u!oy = 

y + 2ti;oi — Swq?/ = (2) 
z + uJqZ = , 

with general solutions [14], [15] 

x{t) =a;o + 2— — 3 ( — — 2yf) ] t — 2— coswoi -I- 2 ( 2— — 3yo ) sincjo* 

y[t) = 2 ( 2yo - — ) + ( 2— - 3j/o ) cos wo< + — sinwoi (3) 
V Wo/ V ^0 / ^0 

z(t) — zq cos Wot H smujot , 

UJQ 

where xo,yo,zo and xo,yo,zo are the initial positions and velocities respectively. 
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Since x(t) contains a term proportional to t, after some time the assumption p Rq 
is no more valid and the above approximation breaks down. However, the divergent 
term can be cancelled by choosing ig = 2a;o2/o- For the LISA case, the constraints of 
rigid and bounded relative motions lead to the solutions [3] 

Xk{t) = ~pa sin [uQt + au] 

Vk (t) = ~ ^/'O cos [uJot + CTfc] 

(t) =-— Pa cos [wo< + CTfc] , 

where (jfc = (1 — fc) 27r/3 and k — 1,2,3 is a label which enumerates the LISA 
spacecrafts. 

2.2. Reseating and expansion of the HCW equations 

For our calculation we rewrite ([2]) by means of the coordinate transformations 

X = x/pii,y = y/pn, z = z/po and i = oj^t. 

x" -2y' ^0 

y" + 2x' ~3y = (5) 



with the notation ' = d/dt. We will refer to ^ as the HCWl equations. We will 
show that the higher order terms in the expansions are negligible for calculating the 
perturbation of the LISA rigid and bounded orbits. 

It is worth noticing that the right-hand side of ^ is zero only at the first order in 
the force expansion. In general, the right-hand side of ([5]) is a polynomial of degree 
n, where n is the order of the expansion |16| . In our case, the order of magnitude of 
the neglected terms is Shcw = Po/^o — 1-9 x 10~^. 



3. Gravitational near field of the EM system 

In our model, the EM system is constituted by 2 point masses: toi ~ 3.0 x 10~^Mq 
(Earth) and m2 ~ 3.7 x 10~^Afo (Moon) located at a constant distance I ~ 
2.57 x 10""^ AU and describing a circular orbit with angular velocity lum = 2tt/Pm, 
where Pm — 29.53 d is the synodic month, around their common barycenter. We 
assume that the barycenter of the EM system makes a circular orbit with radius 1 AU 
around the Sun, i.e. we neglect the eccentricity ~ 0.0167 of the Earth orbit around 
the Sun. 

In addition, we disregard the eccentricity of the Moon orbit around the Earth 
(~ 0.054), its inclination to the ecliptic plane (~ 5.14°), the motion of the perigee 
of the Moon (~ 8.85 yr) and the precession of the Moon orbit plane (~ 18.03 yr) [17], 

m- 

We consider a non rotating reference frame {X, Y, Z} centered on the EM system 
barycenter, with the X axis along the line joining Earth and Moon at t = 0, and the 
Z axis perpendicular to the ecliptic plane (see Figure [2]) . The gravitational potential 
due to the Earth and the Moon at a point {X, Y, Z) is given by 
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U : 
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Zsin(a;A/<) 
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5.-/. Multipole expansion of the EM gravitational potential 

We are interested in the effect of the EM system on the LISA consteUation, located at 
a distance r^. The size I of the EM system is small relatively to r^, l/rg = 7 x 10"'^. 
We therefore expand the total potential in series of em ~ l/^g- 



U{X, Y,Z,t) = y2- el^, Unix, Y, Z, t) , 



(6) 



where C/„ are the well known multipole terms. 
At the zeroth order we have the monopole term 



Uo{X,Y,Z) 



G{mi + 7712) 



(7) 



The first order term, the dipole term, is equal to zero, due to the conservation of the 
linear momentum. The second order term, the quadrupole term, is 



U2{X,Y,Z,t) 



777 1 7772 2 

r„ 



mi + 7772 ^ 



+ y2 _^ ^2 _ 3 [XcOfiLOMt + YsiTLLOMt] 
(A2 + y2 _^ Z2)5/2 



(8) 



and so on. 

Each Un term contains sinusoidal terms as sin kuMt and cos kujMt, with fc = 0, 2, . . . , 77 
for even n and /c = 1, 3, . . . , 77 for odd n. 

Defining com = ojm/^^o — 12.3687 and operating the substitution X ^ x + Xg, Y = 
y + Ug, Z = z in all J7ji, we obtain the EM potential in the HCW frame. The 
corresponding force per unit mass is 

F„ = eEM^Un{x,y,z,t) . 
77! 

The resulting functions F„(x,7/, z,t) are still too complex to be treated analytically, 
but it is worth noticing that the values of the coordinates range within —po and +po, 
while 7'g is about 20 times larger (po/fg = El — 5.5 x lO"^), We can therefore expand 
F„ in terms of x/rg, y/rg, z/rg around the origin of the HCW frame 



E 

F„{x,y,z,t) - V -7 
^-^ 777! 



4 +'^+4' 



The above formula can be expressed in a more useful way (after the rescaling) 
F„(x, ^, - E - E - E FT i'y'i" , 



I ! ^ — ' 7 ! ^ — ' k 

1=0 j=o ■' fc=0 



(9) 



(10) 
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where the an.ijkit) are 

an,tjk{i) = ^ [an,ijk,7nSmmi:uMi + f3n,ijk,mC0srnLdMt] (11) 

with an.ijk,m, f^n,ijk,m numerical coefficients and m = 0,2,...,n for even n and 
m ~ 1,3, ... ,n for odd n. Note that for even n, (|lip contains a constant term (i.e. 
f^n.ijk.o) plus sinusoids, while for odd n it presents only sinusoids {f3n,ijk.i is multiplied 
by cos CoMt). 

Let us estimate the intensity of the acceleration due to the EM system : the most 
important contribution is given by the monopole term {n = 0), and we assume its 
value at the origin of the coordinates as an indicator of its intensity. This term is 

G(mi + mo) q n 

eo = ^ „ , — - ~ 1.31 X 10"-^ ~ 7 X IQ-^encw 
^0 Po 

therefore the EM system influence is about the 7% of the contribution of the Sun 
(Section E2|. 

At the first order in the coordinates the monopole term (normalized to Eq) is 



. Pa{rg -3x1) , - ^PaXgVg 



^9 ^9 

r.2 Q,,2 



Vg I ^ ^PQXgVg . Poirg-Syg) (-^2) 



3 y 3 

rg 



f - Po 
Jz ^ -z — , 

rg 

where i, y, z are the scaled coordinates. 

The most important multipole term is the quadrupole (n = 2) which is mostly 
constituted by the zeroth order term of its expansion in spatial coordinates, or 
equivalently, the force at the origin. This term is periodic with period Pa//2 = tt/w^/, 
and its mean value is 

4 r4 Po mi + TO2 

The monopole force has an identical time-independent polynomial structure as the 
expansion in coordinates of the Sun force. The difference is that this latter starts 
with second order terms [TB], while the EM monopole presents also linear terms and 
a constant term. Moreover, the monopole force is quite smaller than the Sun force. 
The multipole force with even n also contains time-independent terms but the largest 
ones (corresponding to n = 2) have order of magnitude ~ £2 ^ £hcw- 
We can separate the total force in a more useful way: 

(i) A component Fo{x,y, z) independent of time due to the terms of the Sun force 
(oc Ehcw) [13, plus the EM system monopole (oc eo, pO]) with n = 0) plus terms 
contained in the n-even multipole expansion (cx 62, PT|) with m = 0); 

(ii) A component Fi{t) that depends only on time due to the EM system multipole 
terms (cx £2, with i,j,k = and n > 2). Physically this is the force at the 
origin of the HCW system; 

(iii) A component F2(i, y, oc £2 that depends on both time and coordinates (all 
other cases: (fTO|) with n > 2 and m 7^ 0). 
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4. Perturbation of the LISA orbits due to EM system 

If we indicate with roi(i) the unperturbed trajectory of the i*^ LISA test mass 
{i = 1,2,3) and with rii{t) its small perturbation, the difference between two 
perturbed trajectories is simply 

Ar,j(t) EE r,(t) - r^{t) = Aro.y(t) + Ariy(t) , 

where Arojj(i) = roi{t) - roj{t), Ari,^{t) = ru{t) ~ rij{t). 

The perturbation of the relative displacement between the pair i,j of LISA test masses 
can be written as 

where Loij = |Aroy| is the distance between i and j test masses in the imperturbed 
case, and |rii| ^ |roi|. We also define the perturbation to the differential distances 
between each pair of LISA as 

SL,,k{t) = AL„-(t) - ALjkit) (14) 

which represents the variation of Lij{t) — Ljk{t) due to a small perturbation and relates 
directly to the LISA sensitivity curve. 

4-1- Effects of the EM monopole 

At the first order, under the effect of the Sun, the LISA motion is described by (|4]) 
and the EM monopole perturbation is contained in the time-independent force per 
unit mass Fq. 

We write the motion of the fc-spacecraft under the effect of Fq in the following form: 
ffc = fofc + SHCW^ik + £o^2k, where fofc is the unperturbed motion {Qj, rcscalcd), 
and fifc and r2k are the perturbations due to the Sun force terms [16] and the EM 
monopole, respectively. Being sf^cw ^ £o, to calculate f 2fc it is not necessary to know 

fife, [g. 

The equations for r2k = {£2k, y2k, Z2k) are [5], [8] and [9] 

- 2y2fc = fx{xok,yok, zok) 
y'ik + 2i-2fc - 3y2fc = /y(£ofc, yofc, ^ofc) (15) 
z-ik + Z2k = fz{xok,yok, Zok) , 

where {f^, fy,fz) are evaluated along the trajectory fofci using 
The solution can be written as 

X2k = Ci^k + C2,ki+C3MP + (C4,fe + C5,fc^) siu £ + (cg^fe + Cy^fet) COS f 

hk = C8,fc + Cg^ki + (cio,fe + Cii,fcf) siuf + (ci2,fc + Ci3,fc£) COS f (16) 

Z2k = (ci4,fc + ci5,fci) sinf + (cie^fe + Ci7,fcf) cosf , 

where ci^k ■ ■ ■ cn^k are constants that depend on the initial conditions and on the 
geometric parameters {po,Xg,yg). 

The solution contains terms oc t and cx t^, moreover there are also mixed perturbations 
as tsint, i.e. perturbations increases with time. 

In particular, the coefficient of P, cz,k ~ ^•^Xg/{2rg) is positive number and the same 
for all spacecrafts. This means that the entire constellation is "pushed away" by the 
EM system (j-g increases with time). The variation of the LISA arms length can be 
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calculated using Q , ((T5)) and and the result is represented in Figure [3] (right 
panel) . 

The indefinite increasing of the perturbation is not physical because the 

perturbative regime would not be valid anymore, after few years. This is a direct 
consequence of the force linearization: the terms proportional to t sin t and t cos t are 
first-order terms of the real, bounded, solution for small t. Using the F77 code [10], we 
found that the perturbative pulls lead to a complete dismembering of the constellation 
and a successive recombination will occur after several tens of thousand year. In this 
scenario, the distance of each spacecraft from the HCW frame origin ranges from zero 
to 2 AU. Such motion is not a solution of the HCW equations, which are valid only 
if p ^ Rq. An all-time valid solution to our perturbative problem can be obtained 
applying the Lindstcdt-Poincare method [H]. However, as we are interested in the 
LISA motion during a few complete orbits, the Lindstedt-Poincare method is not 
necessary. 

The increase of Vg is shown in Figure |3] (left panel), where its time evolution is 
represented, during the hypothetical first 10 years of the mission. The perturbation 
ALij to the relative motion between the pair i,j of LISA test masses due to the 
monopole perturbation, calculated both analytically (dashed line) and via numerical 
integration [10] (solid line) , is plotted on the right panel of Figure [Sj It is worth 
noticing the good agreement during the 10 years of the LISA mission. 




lime (years) time (years) 

Figure 3: Left panel: Evolution of rg during 10 years: real (solid line) and simplified 
case (Earth describes a circular orbit around the Sun). The periodic component in 
the solid line is due to the eccentricity of the Earth orbit, while the trend is due to 
the EM system influence. i?ight panel: perturbation of the LISA arms length due 
to the EM system monopole effect (rescaled): comparison between the numerical and 
the analytical calculation (solid and dashed lines, respectively). 



4-2. Effects of the EM multipoles 

We now search for a perturbative solution to the HCWl equations in presence of 
n > 2 multipole terms. The intensity of this force is of the order of £2- We have 
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already shown that the muhipole force is composed by a periodic and a polynomial 
component independent of time. The polynomial component is not important, as it 
can be added to the EM system monopole expansion and solved. The solution has 
the same structure as with different coefficients (rescaled of a factor £2/^0) : and 
the motion described in the previous section is therefore a very good approximation 
of the EM system polynomial component influence. 

4-2.1. Fi{t): periodic solutions in nu!i\jt. At the zeroth order {i,j,k = 0) the 
multipole force does not depend on the coordinates. Therefore, for each order n 
we have periodic terms in nujMt. These terms represent the "pure" oscillations of 
a test mass due to the EM system that are not involved with wq harmonics. The 
solutions relative to these frequencies are equal for each test mass, being independent 
of its position (we can interpret this as a common motion). Since we are interested 
in the relative motion of the LISA satellites, we know a priori that these terms are 
subtracted when one measures the distance between two satellites. 
The equations to be solved are 

x" — 2y' = [Oxn SmnLOMt + bxn cos ULOMt] 

n 

[oyn sin nujMt + byn cos nuMt \ (-^T) 

n 

z" + z = [a,„ sin nuiMt + b^n cos ] , 

n 

where axn,yn,zn,bxn,yn.zn are constants. Being nu!i\i 7^ 1 for each n, particular 
solutions can be written as 

X ~ ^ sin nujui + "^o^ nujMt 

n 

y = ^o.'yn sin nLOMt + b'y„ COS ntOMt (Ig) 

n 

z — a'^jj smnuinji + fet„ cos ndiMt 

n 

and the corresponding coefficients are 

^, _ 2bynnLbM + axn(3 + nu)\j) ^, _ - 2a^„ uum + bxn{5 + nLd^j) 



a. 



I _ 2bxn — aynUUlM ,/ _ —2axn — bynnujM , 

I (^zn , / bzn 

^zn ^2 7 zn ^2 7 ' 

nujj^j — 1 ^'^M ^ 1 

Inserting the numerical values, it comes out that the most important contribution is 
due to the 2 cDm frequency and it corresponds to an amplitude of about 1 cm for x, y 
coordinates, while the coefficients a'j,„, 6^„ are all equal to zero. 

In Table [T] we report the coefficients a'xn,ymKn,yn the fluctuations in meters 
(i.e. multiplied by po = 2.89 x 10^ m), relatively to the first six harmonics of the 
fundamental frequency (jJm- 
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Table 1: Rescaled coefficients Po o,xn,yn' PoKn,yn of the sinusoidal terms with 
frequency na)^/ with 1 < ?? < 6. 



n 


Frequency 


Po i't„ 




Po Kn 




po a'yn 




po b'y„ 






[Hz] 




[m] 






[m] 






[m] 




[m] 




1 


3.920> 




-4.5; 


<io- 


-5 


-1.8x 


:10- 


-4 


6.5x10- 


-5 


-4.4x10- 


-5 


2 


7.840 > 




4.4; 


<io-^ 


-3 


8.5x 


:10- 


-3 


-5.7x10- 


-3 


3.9x10- 


-3 


3 


1.176> 


(10-" 


-2.1; 


<io- 


■d 


-2.7x 


:10' 




1.9x10- 




-1.8x10" 


-5 


4 


1.568> 


<io-*' 


1.1; 


<io- 


-7 


l.Ox 


:10- 


-7 


-7.6x10- 


-8 


9.8x10- 


-8 


5 


1.960> 


(lO-** 


-6.4; 


<io- 


-10 


-3.9x 


:10- 


-10 


3.0x10- 


-10 


-5.7x10- 


-10 


6 


2.352> 


(lO-** 


3.8; 


<io- 


-12 


1.5x 


:10" 


-12 


-1.1x10- 


-12 


3.4x10- 


-12 



4-2.2. F2{x, y, z,t): periodic solutions in {nuj^ ±m)t. In Section 14.2.11 we showed 
that the Fi term corresponds to the particular solutions of (|17p and that Fi does not 
affect ALij being independent of coordinates. 

The motion associated with the coordinate-dependent term F2, will be different 
between each pair of test masses, and so the relative displacements ALij will be 
different from zero. 

The most direct approach to solve the equation of motion is to write HCWl equations 
with F2 evaluated along the unperturbed trajectories given by Q. Thus, we have 
obtained only the amplitudes relative to frequencies (nujM i 1), which represent the 
main effect of F2. In order to have the complete spectrum (jiujm i to) we should 
consider the solutions of HCW with the complete expansion of the Sun force per unit 
mass (Section 12. 2|) . 

The solution can be written as sum of sinusoidal terms with frequencies {nujM ±1). In 
addition, using (jl3p . a similar relation can be written also for the perturbation ALy 
that in SI units reads 

1 

AL.y(i) = poX! X! [a„„i,ij sinw„„ii + 6„^,.y cosw„m,i] , (20) 

n m— — 1 

where w„m = nujM + 2muQ, with n = 1,2,..., -l-oo and to = —1, 0, -1-1. 

The quantity SLiji^, which can be directly compared with the LISA sensitivity curve, 

reads 

1 

5Lijk{t) = Po [anm.ijkSVI^i^nmt + bnm.ijk cos LOnmt] , (21) 

n m— — 1 

where O^nyn.ijk — ^nva.ij ^nmjk and &7ijn,zjA: — ^nm^ij ^nni.jk- 

The coefficients of (PUI) and (PT|) arc reported in Table [5] and Table [31 respectively, 
relatively to a certain number of frequencies. 

In Figure |3] wc plot the SL123 amplitudes (filled circles) superimposed to two LISA 
sensitivity curves (straight lines) corresponding to integration times of 12 days (the 
upper one), period below which there should not be disturbances or if present should 
be removable and 1 year (the lower one), respectively. The straight lines were obtained 
by extrapolating the LISA sensitivity curve down to 10~^ Hz, as discussed in [I]. 
Each amplitude is subject to time variations, because i) the Earth and the Moon orbits 
are not circular; ii) the orbital plane of the Moon is slightly inclined; and iii) the LISA 
constellation is not rigid. 
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We estimated uncertainties for Xg,yg,l and po by taking into account all these effects; 
in particular, we assumed the error relative to po equal to the amplitude of the Sun 
induced breathing, i.e. ~ 2% of po E] • 

The resulting relative uncertainties of the SLijk amplitudes are ~ 30%. 

0.01 I 1 1 1 1 — I — I — I 1 1 

= S = 

0.001 r * s 

0.0001 1-j -1 

E 1 E 



10-5 ^ I 
E i 




0.4 0.6 0.8 1 2 
frequency (10-^ Hz) 



Figure 4: The filled circles give the amplitudes of the differential relative displacement 
(5Li23 between the LISA test masses 1,2 and 2,3 induced by multipoles of the EM 
system; relative errors are also plotted. The two straight lines represent the LISA 
sensitivity curve for 12 days (upper straight line) and 1 year (lower straight line) of 
integration time. 



5. Effect of Venus and Jupiter on LISA motion 

To make a comparison with the effect on the LISA motion due to EM system effect, 
we have numerically evaluated the effect of Venus, Jupiter and the EM system (its 
monopole contribution) using a F77 code implementing the algorithm described in 
[10]. In practice, we independently computed the LISA motion under the effect of 
Sun+ Venus, Sun+ Jupiter and Sun+EM system and we subtracted, to each of them, 
the unperturbed motion due to the Sun. In Figure [5] we plot the SL123 perturbations 
(in km) once the modulation due to the Sun has been subtracted. We have found that 
the monopole contribution of the EM system (solid line) is much larger than that of 
Venus (dotted line) and Jupiter (dashed line). Figure [5] also shows the onset of the 
resonant effect of the EM system monopole after 2-3 years. 
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Table 2: Coefficients relative to ALij{t). Amplitudes and frequencies are in SI units 
and the overall accuracy on amplitudes is ~ 30%. 



n, m) 


Frequency 
[Hz] 




[m] 




^Tim,12 

[m] 




finm,23 

[ni] 




^7H71,23 

[m] 




71 771, 31 

[m] 




^71777 , 31 

[m] 




L-1) 


3.286x10" 


7 


5.0x10" 


-.5 


-8.0x10" 


-6 


-1.8x10" 


-5 


4.8x10" 


-r. 


-3.2x10" 


-5 


-4.0x10" 


-5 


1, 0) 


3.920x10" 


7 


6.8x10" 


-6 


1.6x10" 


-5 


6.8x10" 


-6 


1.6x10" 


-r. 


6.8x10" 


-6 


1.6x10" 


-5 


1, 1) 


4.553x10" 


7 


1.9x10" 


-5 


-1.4x10" 


-5 


2.7x10" 


-6 


2.4x10" 


-5 


-2.2x10" 


-5 


-9.5x10" 


-6 


2rl) 


7.206x10" 


7 


-2.1x10" 


-3 


-2.1x10" 


-4 


1.2x10" 


-3 


-1.7x10" 


-3 


8.8x10" 


-4 


1.9x10" 


-3 


2, 0) 


7.839x10" 


7 


-6.4x10" 


-4 


-7.9x10" 


-4 


-6.4x10" 


-4 


-7.9x10" 


-4 


-6.4x10" 


-4 


-7.9x10" 


-4 


2, 1) 


8.473x10" 


7 


-3.9x10" 


-4 


5.0x10" 


-4 


-2.3x10" 


-4 


-5.9x10" 


-4 


6.2x10" 


-4 


9.1x10" 




3,-1) 


1.113x10" 


6 


8.9x10" 


-6 


2.6x10" 


-6 


-6.7x10" 


-6 


6.4x10" 


-6 


-2.2x10" 


-6 


-9.0x10" 


-6 


3, 0) 


1.176x10" 


6 


3.6x10" 


-6 


3.0x10" 


-6 


3.6x10" 


-6 


3.0x10" 


-6 


3.6x10" 


-6 


3.0x10" 


-6 


3, 1) 


1.239x10" 


6 


1.1x10" 


-6 


-2.2x10" 


-6 


1.4x10" 


-6 


2.0x10" 


-6 


-2.5x10" 


-6 


1.6x10" 


-7 


4,-1) 


1.505x10" 


6 


-4.4x10" 


-8 


-2.2x10" 


-8 


4.2x10" 


-8 


-2.7x10" 


-8 


2.6x10" 


-9 


5.0x10" 


-8 


4. 0) 


1.568x10" 


6 


-2.3x10" 


-8 


-1.3x10" 


-8 


-2.3x10" 


-8 


-1.3x10" 


-8 


-2.3x10" 


-8 


-1.3x10" 


-8 


4, 1) 


1.631x10" 


6 


-3.0x10" 


-9 


1.2x10" 


-8 


-8.5x10" 


-9 


-8.4x10" 


-9 


1.2x10" 


-8 


-3.2x10" 


-9 


5,-1) 


1.896x10" 


6 


2.3x10" 


-10 


1.8x10" 


-10 


-2.7x10" 


-10 


1.1x10" 


-10 


3.7x10" 


-11 


-2.9x10" 


-10 


5. 0) 


1.960x10" 


6 


1.5x10" 


-10 


5.4x10" 


-11 


1.5x10" 


-10 


5.4x10" 


-11 


1.5x10" 


-10 


5.4x10" 


-11 


5, 1) 


2.023x10" 


6 


3.9x10" 


-12 


-6.5x10" 


-11 


5.4x10" 


-11 


3.6x10" 


-11 


-5.8x10" 


-11 


2.9x10" 


-11 


6,-1) 


2.288x10" 


6 


-1.2x10" 


-12 


-1.3x10" 


-12 


1.7x10" 


-12 


-3.8x10" 


-13 


-5.4x10" 


-13 


1.7x10" 


-12 


6, 0) 


2.320x10" 


6 


-9.9x10" 


-13 


-1.7x10" 


-13 


-9.9x10" 


-13 


-1.7x10" 


-13 


-9.9x10" 


-13 


-1.7x10" 


-13 


6, 1) 


2.415x10" 


6 


5.0x10" 


-14 


3.7x10" 


-13 


-3.5x10" 


-13 


-1.4x10" 


-13 


3.0x10" 


-13 


-2.3x10" 


-13 



Table 3: Coefficients relative to SLijk{t). Amplitudes and frequencies are in SI units 
and the overall accuracy on amplitudes is ~ 30%. 
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Figure 5: 5Li2Z perturbation (in km) due to EM system (solid line), Venus (dotted line) and Jupiter 
(dashed line). &L123 perturbation (in km) due to monopole contribution of the EM system (solid 
line), Venus (dotted line) and Jupiter (dashed line). See Section [5] 



6. Conclusions 

We calculated the relative motion of LISA test masses due to the effect of the EM 
system monopole and we found AL^j ~ 3 x 10^ km after a period of 10 years. We 
also found that the Jupiter and Venus influences are at least 10 times smaller than 
the EM system one. 

The perturbations of the differential relative motion of LISA test masses dLijk are in 
the 10~^ to 10~^ Hz decade. This is a very low frequency range in which LISA residual 
acceleration noise may be much larger than extrapolated on grounds of known effects 
in [T]. However, it might not be completely hopeless to get an interesting sensitivity 
also at such low frequencies, provided that one will have to face the problem of i) 
loss of signal coherence over a time scale of one month; and ii) low frequency range 
calibration. 

As discussed in Pollack it is possible to extract a signal from the LISA data, 
even in presence of disturbances. These latter arise due to environmental effects, 
such as cosmic rays induced by solar flares, and the telecommunication antenna which 
periodically has to be rotated. Pollack showed how these disturbances can be identified 
and subsequently removed from the data even at low frequencies. Assuming, for 
instance, that a disturbance appears every 19 days (see Table 6 in [20]) the resulting 
error in the signal frequency of 3 x 10~^ Hz is only ~ 1.3 x 10~^ Hz. Thus it should not 
be a problem to extend to such low frequencies the calibration from the verification 
binaries, and, by using their signals, ensure the continuity of data over time spans 
of many weeks. Still the signal from the EM system, as understood at the level of 
accuracy given in the present paper, can give a relevant additional crosscheck to such 
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an extension of the calibration. It thus may help in improving our knowledge of the 
LISA acceleration noise at very low frequencies and contribute to extend to such low 
frequencies the capabilities of LISA. 
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